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Anisotropic order parameter 
in Li-intercalated layered superconductor Lia,ZrNCl 
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In this study, it is shown that in a layered nitride superconductor, i.e., Li^^ZrNCl (0.07 < a; < 0.21), 
the superconducting order parameter is highly anisotropic in a sample with x — 0.12, as inferred 
from both the temperature and the magnetic field dependences of the muon depolarization rate 
((Ts, proportional to the superfluid density). Moreover, the tendency of strong anisotropy with an 
increase in x is indicated by the T dependence of as. These observations are in good agreement with 
the recent theory that predicts the development of anisotropy in a d + id' gap upon carrier filling 
to the bands with disconnected Fermi surfaces on a honeycomb lattice. 

PACS numbers: 74.70.-b, 74.20.Rp, 76.75. +1 



Layered nitrides such as /3-MNCl (with M = Zr, 
Hf) are attracting considerable attention since they ex- 
hibit superconductivity upon intercalating alkaline met- 
als [U, 01 . While these nitrides have relatively high su- 
perconducting transition temperatures [Tc ~ 15 (25) K 
with M — Zr (Hf)], the density of states at the Fermi 
level [A^(0)] is reported to be considerably lower than 
that of other superconductors having similar Tc, as in- 
ferred from the results of the specific heat measurement of 
Lio.iaZrNCl i and magnetic susceptibility measurement 
of Lio.48(THF)yHfNCl (where THF refers to tetrahydro- 
furan) [J]. The small A''(0) is in excellent agreement with 
the prediction of the theoretical investigation^. Such 
situation naturally leads to the question of whether or 
not the superconductivity in layered nitrides is fully un- 
derstood on the basis of the conventional BCS theory 
with electron-phonon coupling. Moreover, it has been 
shown that one of these nitrides, i.e., Li^jZrNCl, exhibits 
further anomalies that are not expected to be present in 
the simplest situation presumed by the BCS theory. 

Lia;ZrNCl has a lamellar structure consisting of alter- 
nating stacks of Zr-N double honeycomb layers and insu- 
lating CI bilayers. Li atoms are intercalated into the van 
der Waals gap of the CI bilayers to supply electrons to 
the conducting Zr-N layers. In addition to the above 
mentioned small N{0), it has been also revealed that 
the electron-phonon interaction is too weak to explain its 
high TcIIjSS]- general, in the case of two-dimensional 
electronic systems such as /3— ZrNCl, A^(0) may be only 
weakly dependent on band filling. Therefore, provided 
that superconductivity is explained by the conventional 
BCS theory, Tc would not vary with the Li concentration 
(x). However, the fact is that while Tc is independent of 
doping for x > 0.15, it increases steeply below x ~ 0.12, 
reaching a maximum (Tc = 15.2 K at a; = 0.06) and 
then suddenly transitioning into an insulating state for 
X < 0.05 0. It might be worth noting that this tendency 



of Tc to be high at a low carrier density x is opposite to 
that of underdoped cuprates. 

Another anomaly is reported in the magnetic field de- 
pendence of an electronic specific heat (Sommerfeld) co- 
efficient 7 in the mixed state of Lip. 12 ZrNCl. While 7 is 
expected to be approximately proportional to the num- 
ber of flux lines, and accordingly to the magnetic field in 
conventional BCS superconductors (i.e., 7 ~ 7„i?/i?c2 oc 
H/^o, where 7„ is the electronic specific heat in the 
normal state. H is the magnetic field, Hc2 is the upper 
critical field, and (f>o = 2.07 x 10~^^ T-m^ is the quan- 
tum flux), it increases with a gradient much steeper than 
jnlHc2 in Lio.i2ZrNCl, approaching slI H ^ 0.4i?c2 
Q . This strongly suggests the occurrence of field- induced 
quasiparticle excitation that is not expected for the su- 
perconducting order parameter described by isotropic s- 
wave paring with a unique gap energy. 

In this paper, we describe our muon spin rotation 
(/iSR) study on Li^; ZrNCl over a range of Li content near 
metal-to-insulator transition (0.07 < a; < 0.21). We show 
that in a sample with x = 0.12, the temperature and 
magnetic field dependences of the muon spin depolariza- 
tion rate [as = Us{T,H), which is proportional to the 
superfluid density n^] is in complete agreement with the 
anomalies observed in the bulk properties, providing mi- 
croscopic evidence for the anisotropic order parameter. 
Furthermore, a clear tendency of increasing anisotropy 
as a function of x is inferred from the magnitude of the 
gap ratio (2A/fcBTc). These observations support the re- 
cent theory that predicts the occurrence of anisotropic 
d -\- id! pairing in Li^; ZrNCl and the effect of electronic 
correlations on the Fermi surface specific to Li^: ZrNCl, 
which develops with band filling 0] . 

Conventional ^SR measurements were performed on 
M15 and M20 beamfines of TRIUMF, Canada. A ^SR 
apparatus with high time resolution was used to mea- 
sure the time-dependent positron decay asymmetry un- 
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FIG. 1: (Color online) (a) X-ray diffraction spectra for the 
present samples obtained at SPring-8 (X-ray energy of 15.6 
keV). The [104] line observed in the pristine compound is 
absent in Li-doped samples, (b) Lattice parameters (open 
circles show present data, and filled points are after Ref. 

der a transverse field (TF) up to 2 T. LiajZrNCl samples, 
alined along the c axis, were loaded on a He gas-flow 
cryostat, and they were field cooled to a target temper- 
ature to minimize the effect of flux pinning. Details of 
sample preparation are described in the earlier report, 
where special precaution was taken to maintain the ho- 
mogeneity of samples 0. The superconducting volume 
fraction as well as [= 15.1 K {x = 0.07), 14.0 K 
{x = 0.08), 12.6 K (a; = 0.10), 12.5 K (a; = 0.12) and 
11.6 K (x = 0.21)] determined from magnetization mea- 
surements were in excellent agreement with the previous 
result . The concentration of lithium was determined 
from inductively coupled plasma spectroscopy. As shown 
in Fig. [TJ the homogeneity of samples was confirmed by 
carrying out high-resolution powder X-ray diffraction at 
SPring-8. The mean free path I of these samples were 
estimated to be 10-16 nm (increasing with a;)[ly|, which 
is approximately equal to the coherence length (^ = 8-13 
nm). Although this would not be in the clean limit, the 
condition I > ^ is always satisfied, and thus, the situation 
is consistent with the observed influence of anisotropic 
gap on CTs. 

^SR is an effective microscopic technique for measur- 
ing the magnetic penetration depth (A) in type II super- 
conductors. It is reasonably presumed that implanted 
muons are distributed randomly over the length scale of 
flux line lattice (FLL), probing local magnetic fields at 
their respective positions. Then the transverse muon spin 
precession signal consists of a random sampling of inter- 
nal field distribution B{r), such that 

/oo 
cos(7^Bt + <f>)n{B)dB, 
'OO 

n{B) = {d[B{r) - B])r, 

where 7^ = 27r x 135.53 MHz/T is the muon gyromag- 
netic ratio, n{B) is the spectral density for the internal 
field defined as a spatial average (()r) of the delta func- 
tion, and (j) is the initial phase of rotation. Hence, n{B) 
can be obtained from the real amplitude of the Fourier 
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FIG. 2: (Color online) Examples of TF-/iSR spectra with 
X = 0.21 (a) and 0.12 (b) shown in a rotating reference frame 
frequency of 20 MHz and 40 MHz, respectively. Circles repre- 
sent the spectra above Tc, while triangles are those obtained 
below Tc. Solid lines are fits obtained using Eq. ([T|. 

transform of the TF-/xSR signal. In the case of a rela- 
tively long magnetic penetration depth (A > 300 nm), 
the Gaussian distribution is a good approximation for 
n{B), yielding 

P^{t) ~ exp(-CTf tV2) cos(cjot + (I)), 

where CTs is obtained from a second moment of the field 
distribution (= 7^-^/ {[B{r) — S]^}), and wq = 7^ Bo with 
Bq ~ ^J■oH. Here, provided that the clean limit is 
achieved, A is related to the superconducting carrier den- 
sity Hg as follows: 

1 n^e^ 

Fig. [2] shows some examples of the TF-/.tSR time spec- 
tra observed in the samples with x = 0.21 and 0.12 un- 
der a field of 0.15 T and 0.3 T, respectively (shown in a 
rotating- reference- frame for visibility) . While both spec- 
tra show a slow Gaussian damping above Tc because of 
random local fields from nuclear magnetic moments, fur- 
ther depolarization resulting from the formation of a flux 
line lattice is observed upon superconducting transition. 
A clear trend toward a high depolarization rate is ob- 
served for a high Li concentration x, which indicates an 
increase in Ug with x. We also note that an additional 
component that exhibits fast depolarization (^2 MHz) is 
observed in samples with x = 0.08, 0.10, and 0.12. This 
observation is tentatively attributed to muonium (Mu) 
that may be formed when muons are stopped near the 
insulating CI bilayers. Considering a background con- 
tribution from muons stopped in the sample holder, we 
used the following function for the analysis of the /iSR 
time spectra by curve fitting; 

AoP^it) = cxp{-<jlt^)[AgP^{t) + Afe~^' cos(cjot + 4>)] 
+ A), exp(-alt'^) cos {ujbt + (/)), (1) 
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where Aq is the total positron decay asymmetry (^--^0.2), 
CTn is the depolarization rate attributed to nuclear mag- 
netic moments, As is the partial asymmetry of a super- 
conducting fraction, A{ is that of the component related 
with Mu formation showing depolarization at a rate A, 
Ah (= ^0 — — < 0.01) and Uh are the amplitude 
and central frequency of the background. It was found 
that A{ and A were mostly independent of temperature 
and magnetic field, and their fractional yield {Af/Ao) 
was 0.26(1), 0.169(1), and 0.163(1) for x = 0.08, 0.10, 
and 0.12, respectively. The influence of flux pinning on 
ffs was confirmed to be negligible for x = 0.10 and 0.21, 
as inferred from the nearly field-independent Us at low 
magnetic induction (< 0.2 T). This also confirms the 
formation of three-dimensional vortices as expected from 
relatively large coherence length (^ = 8-13 nm) over the 
Cl-bilayer thickness (^1 nm). 

Figure [3] shows a plot of Cs against temperature for 
the samples in which CTs has been deduced with suffi- 
cient precision; it was found that rig was too small for 
samples with x < 0.10. Solid curves are the best fits 
obtained by applying the conventional BCS theory for 
s-wave symmetry with a single-gap (s-BCS theory) pT| . 
Here, assuming that electron-phonon coupling plays a mi- 
nor role in Li^^ZrNCl, we use the s-BCS theory to deter- 
mine the curve fit as an effective model extended to the 
case of anisotropic order parameters (having dips, nodes, 
or a secondary gap in multiple bands). This extension 
is achieved by allowing the gap ratio (2A/fcBrc) to vary 
freely, where the energy gap A is regarded as a mean 
value. A, averaged over the momentum space and rele- 
vant bands (while is determined by the maximum of 
A). Then, we find that the deduced gap ratio (2A/A:bTc, 
shown in the inset of Fig. [3)) decreases with an increase 
in x. Assuming that decreases with an increase in 
X, we find that the decrease in A is steeper than that 
in decrease in Tc- While this tendency is not easily un- 
derstood in the framework of electron-phonon coupling, 
it is understood in the extended model as an indication 
toward strong anisotropy with high Li content. 

While we cannot distinguish the origin of anisotropy 
between one-band and multiband scenarios solely from 
the behavior of supcrfluid density, the presence of a sec- 
ondary energy scale in the energy gap is inferred from 
the analysis using a phenomenological double-gap model 
for s-wave symmetry [l2j, such that 

a^iT) ^ [1 - wSaiAi,T) - (1 ~ w)Sa{A2,T)] , 
2 



6a{S, T) 



/(e,r)-[l-/(e,r)]&. 



/(e, T) = [l + exp (.J^TWY/kBT 



where f{e,T) is the Fermi distribution function, A(r) 
is the standard BCS gap energy, and w is the fractional 
weight of the i = 1 component. Although it is not easy 
to determine all the parameters using the curve fits of 
data shown in Fig. |31 A2 for x = 0.12 can be determined 
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FIG. 3: (Color online) Temperature dependence of CTs for sam- 
ples with X = 0.10, 0.12, 0.21. Solid cmves are fits obtained 
using the s-BCS theory as an effective model. Longitudinal 
and horizontal axes are normalized by crs(O), Tc for compar- 
ison. Inset shows the x dependence of a fitting parameter 
2A/A;bT'c, with A interpreted as a mean value. 



on the basis of an assumption that Ai corresponds to 
A(T = 0) determined by a jump in the specific heat near 
Tc (2.64 meV) yielding excellent fit with the data 
obtained when 2A2/kBTc = 2.4(4) [A2 = 1.3(2) meV, or 
A1/A2 = 2.0(3)] and w = 0.59(9). 

The presence of a secondary energy scale in the or- 
der parameter is further suggested by the magnetic field 
dependence of (Tg measured at 2 K for the same sample 
in which an anomalous behavior of 7 has been reported 
3 . As shown in Fig. HI it is observed that in the limit 
of H/Hc2{= ^) — > 1, CTg exhibits a trend of asymptotic 
conversion to a value (ctc). In the case of muonium for- 
mation, we attribute CTc to an artifact resulting from an 
incomplete separation of As and Af in the curve fit that 
might have led to residual depolarization, and we model 
the field dependence as follows: 



(7s{h) = Z(7Q{hl) + (1 - Z)(7a{h2) + (Tc 

ao{hi) = 0.0274 x ^^(1 - h,) [l + 3.9(1 - h,f] 



(2) 

1 

(3) 



where hi = H/H^2 ^ is the secondary upper critical 
field (a parameter corresponding to A2), z is the rela- 
tive weightf of (Ts(^i), and Eq. Q is an approximated 
expression for the field dependence of erg for a single -gap 
case (with the Ginzburg-Landau parameter k ^ 1) [l3|. 

In the curve fit, h'^J^ was fixed to the reported value 
(5 T, determined by the specific heat 3). This model 
reproduces our data excellently, yielding z = 0.22(3) and 
1.2(1) T. A good agreement between the ratios 

,^2^ = 4.2(2) and {A1/A2Y = 4.1(5) is perfectly 
in line with the general relation Hc2 oc A^. (Note that 
the ratio A1/A2 is also in excellent agreement with that 
evaluated from the T dependence.) More interestingly, 
the value of Hc2 corresponds to the field toward which 
7 exhibits a steep increase and then gradually saturates 
to 7„ for H^l^ < H < h'^J^ 0. This, together with 
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FIG. 4: (Color online) Field dependence of CTs in Lia:ZrNCl 
with X = 0.12. Solid curve represents the best fit obtained 
by Eq. ((Sjl, and the dashed curve corresponds to the case of 
z = 1. Inset shows magnetic penetration depth deduced from 
the same analysis. 



the T dependence of a^, strongly suggests that the su- 
perconducting order parameter is highly anisotropic in 
the sample with x = 0.12 that is characterized by a sec- 
ondary energy scale A2 ~ 0.5Ai, where A2 may result 
from dips in a single-band gap or double-gap structure. 
Here, it would be worth quoting a prediction of band- 
structure calculation that carrier filling to the secondary 
band would not occur until x exceeds ~0.3, where a large 
increase in N{0) is expected d, [l!]- Absence of such an 
increase in the recent specific heat measurement over the 
relevant range of x supports the one-band scenario [l^ , 
suggesting that A2 corresponds to the minimal gap en- 
ergy (at the dips) in the order parameter. 

According to a recent theory based on the Hubbard 
model considering disconnected Fermi surfaces on a hon- 
eycomb lattice [8|, spin fluctuation enhances Tc over a 
low X region. Since spin fluctuation-mediated supercon- 
ductivity needs to have sign change in the gap func- 
tion, the structure of the relevant Fermi surface leads 
to the prediction of d + id' pairing symmetry as the 
most probable candidate. The decrease in Tc with an 
increase in doping is explained by an increase in three 
dimensionality and associated reduction in the relative 
volume in the Brillouin zone, where the pairing interac- 
tion is strong (corresponding to wave vectors that bridge 
the opposite sides of each pieces of the Fermi surface). 
Moreover, the anisotropy attributed to the dips in the 
order parameter develops for x > 0.11 and that the min- 
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imum of gap energy is reduced to Ami 
X = 0.16 [9|. This behavior is in qualitative agreement 
with the tendency suggested by the x dependence of the 
gap ratio 2A/kBTc deduced from the T dependence of as, 
where the secondary energy scale suggested in the case of 
X = 0.12 (A2 ~ O.5A1) may correspond to the minimal 
gap (Aniin A2). Then, it is likely that the decrease in 
the gap ratio with an increase in x reflects a decrease in 
Amin in the d + id' gap. 

Finally, we discuss the behavior of (Ts{T 0) as a 
function of Tc. As shown in Fig. [SJa), Tc increases with 



a decrease in as- It is inferred from our data shown in 
Fig. ^h) that ag is mostly proportional to x. Thus, the 
CTs (x) dependence of Tc is in excellent agreement with the 
earlier result that is shown in Fig.[5l^c) [7[, which strongly 
supports the high quality of the investigated samples. 
Meanwhile, this is in marked contrast with the as depen- 
dence of Tc reported in an earlier literature [l6j . where 
the authors maintain a linear relation common to that 
observed in underdoped cuprate superconductors (indi- 
cated by shaded area) based on their result obtained for 
X = 0.17 and 0.4 (open triangles). In this regard, we 
point out the fact that the Li^jZrNCl samples studied in 
Ref. [1^ do not follow the Tc versus x relation observed 
in our samples. In particular, their sample with x = 0.17 
seems to exhibit relatively high Tc corresponding to that 
oi X — 0.07-0.08 in our sample. Considering the diffi- 
culty associated with obtaining a uniform specimen for a 
low Li concentration Q, one might speculate that their 
sample (particularly for x = 0.17) might have had prob- 
lems concerning homogeneity. The reported small ^SR 
asymmetry (:^0.1 [l^) might be further evidence for this 
speculation. In any case, our result shows that the rig de- 
pendence of Tc is markedly different from that observed 
in the case of underdoped high Tc cuprates. 
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FIG. 5: (Color online) (a) Tc as a function of as(T 0). 
Solid circles show the present data, and open triangles are 
quoted from ^ (only for Li^ZrNCl with x = 0.17 and 0.4, 
where the latter corresponds to the point with lower Tc). The 
hatched area represents the empirical linear relation found in 
underdoped cuprates [31 • (b) x dependence of cts (the present 
data) with the dashed line representing a linear relation. This 
line crosses zero at 2; ~ 0.05, where metal-insulator transition 
occurs, (c) X dependence of Tc (after Ref.Q]). 

In conclusion, it is microscopically shown by ^SR mea- 
surements that the superconducting order parameter in 
LiajZrNCl is strongly anisotropic in a sample with x = 
0.12, as inferred from the temperature and field depen- 
dence of <Ts- Therefore, the origin of anomalous behav- 
ior in the Sommerfeld coefficient observed in the specific 
heat measurements is attributed to the strong anisotropy 
in the order parameter that might be characterized by 
a secondary energy scale A2 (where A2 ~ 0.5Ai). The 
temperature dependence of as over a range of x from 0.10 
to 0.21 implies a tendency for the occurrence of a weak 
pairing interaction and strong anisotropy in the order pa- 
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rameter with an increase in x. The latter feature as weU 
as the relatively high at a low Li concentration sup- 
ports the important role of electronic correlation for the 
mechanism of superconductivity in Lij;ZrNCl predicted 
by the recent theory. 
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